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Resistivity is an intrinsic feature that specifies the electrical properties of a material and depends on 
electron-phonon scattering near room temperature. Reducing the resistivity of a metal to its potentially 
lowest value requires eliminating grain boundaries and impurities, but to date few studies have focused on 
reducing the intrinsic resistivity of a pure metal itself. We could reduce the intrinsic resistivity of 
single-crystal Ag, which has an almost perfect structure, by impurity doping it with Cu. This paper presents 
our results: resistivity was reduced to 1.35 ufi-cm at room temperature after 3 mol% Cu-doping of 
single-crystal Ag. Various mechanisms were examined in an attempt to explain the abnormal behavior. 

Metals are considered to be ideal materials for electronics, because of their superior electrical properties 1 . 
Silver (Ag) is the best conductor at room temperature, with a resistivity of 1.59 ufl-cm at 293 K, which 
maybe expressed as 108% according to the International Annealed Copper Standard (IACS). The high 
conductivity of metals, especially Ag, makes them ideal for many industrial power applications, in which excellent 
printing, thermal, electrical, wire-bonding, and soldering properties are a prerequisite". Many researchers have 
recently tried to develop advanced metal materials with greater power consumption and higher efficiency for the 
aerospace and medical industries 714 . Although superconductors are ideal for this purpose, their low critical 
transition temperature (T c ) remains a limiting factor for practical applications. Thus, considerable research 
has focused on developing new materials with higher conductivity at room temperature (RT) 15 " 17 . 

The resistivity of a material increases as a result of imperfections, such as defects, impurities, grain boundaries, 
and dislocations 18 . The resistivity of single-crystal copper (Cu) (1.52 ufl-cm at 293 K, 1 13.4% IACS) is 12% lower 
than that of polycrystalline Cu, which is even lower than that of polycrystalline Ag, due to the elimination of grain 
boundaries via the crystal growth process 19 . Another approach to reducing resistivity involves the simultaneous 
application of heat and pressure after crystal growth, which can enhance the crystallinity of the metal and thereby 
significantly reduce resistivity (e.g., for Cu, a 14% reduction in resistivity: 1.472 ufl-cmat293 K, 117.1% IACS) 20 . 

The resistivity of alloys is generally higher than that of pure metals, and increases with impurity contamina- 
tion 21,22 . However, we report here a surprising decrease in resistivity by addition of a certain amount of Cu 
impurity to single-crystal Ag, which is an ideal conductor at room temperature. Copper wire of high purity has a 
conductivity of slightly over 103% IACS, which is close to the value expected for Cu with no impurities. Here, we 
report an uncharacteristic increase in conductivity, to 127% IACS, as a result of Cu impurities in single-crystal Ag. 

Single-crystal Cu, single-crystal Ag, and Ag!„ x Cu x mixed crystals (x = 1, 2, 3, 5, and 10 mol%) were grown 
using the Czochralski method (all equipment was supported by Crystal Bank at Pusan National University). 
Figure 1 shows an example of an Ag 0 95 Cu 0 05 mixed crystal. These pole-figure images, obtained using X-ray 
diffraction (XRD) (Fig. lb, c), confirm the crystallinity of the Ag!_ x Cu x crystals (additional pole-figure images of 
other Ag!_ x Cu x mixed crystals are provided in Supplementary Information Fig. 1). The crystal quality of the Ag- 
Cu mixed crystals was slightly lower than that of pure single-crystal Cu and Ag, but the mixed crystals exhibited 
well-defined growth along the (111) direction, up to a Cu-alloy ratio of 10 mol%, and most of the copper atoms 
were substituted into Ag sites. Samples for resistivity measurements were prepared by wire cutting using electrical 
discharge machining (EDM), as shown in the inset to Figure le. 
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Figure 1 | Grown Ag-Cu mixed crystal and its fabrication. 

(a): Photograph of the 5% Cu-doped Ag mixed crystal, grown using the 
Czochralski method, (b) and (c): 2.5D and 2D pole-figure images obtained 
from the metal mixed crystal, (d) and (e): Sample-cutting image using a 
wire electrical discharge machine (EDM). Inset to (e): Sample for electrical 
resistivity measurements using the four-point probe method. All 
photographs were taken by S.Y.J, and the co-authors. 

Compositional analysis was carried out using glow- discharge 
spectrometry (GDS). We confirmed that the Cu atoms were homo- 
geneously distributed with the intended amount. Concentrations of 
impurities other than Cu were negligible, compared with that of Cu. 
The distribution of unintentional impurities, other than Ag and Cu, 
are shown in the Supplementary Information section (Fig. S2 and 
Table 1). 

Figure 2a shows the electrical resistivity for various Ag/Cu materi- 
als at room temperature. The resistivity of metals is known to 
increase with the addition of impurities, and is also higher in mixed 
crystals or alloys than in pure crystalline materials. 

However, we observed, unexpectedly, that resistivity decreased 
after doping with small amounts of Cu, up to 3 mol%. The lowest 
resistivity was obtained in the 3 mol% Cu (^)-doped Ag mixed 
crystal; its resistivity was 1.35 ufl-cm, which is 9.4% lower than that 

of single-crystal Ag (#) (1.49 ufi- cm). The resistivities near 3 mol% 
(±0.1, 0.2%) of Cu were slightly higher than that of 3 mol% (inset, 
Fig. 2(a)). The resistivity of greater than 5 mol% Cu-doped Ag was 
slightly higher than that of single-crystal Ag. Resistivity increased 
slightly at greater than 5 mol% Cu-doped Ag. The resistivities of the 
5 mol% Cu (M)- and 10 mol% Cu ( )-doped Ag crystals were 1.51 
and 1.66 uQ-cm, respectively, slightly higher than that of pure sin- 
gle-crystal Ag. This is the first report of this uncharacteristic drop in 
resistivity due to Cu doping. 

In contrast, the alloy (not grown as a single crystal) exhibited an 
increase in resistivity as a result of Cu addition. The alloy was 



prepared by furnace cooling from the melt state. The 3 mol% Cu- 
doped alloy (O) had a resistivity of 1.76 |ifl-cm, a 23.3% increase 
over that of the Ago.97Cuo.03 mixed crystal. 

The different consequences of Cu doping in the alloy and the 
single crystal may be attributable to the short-range order of Cu in 
the mixed crystal and the reduced number of grain boundaries. The 
liquid cooling rate is much slower for crystal growth, resulting in an 
increase in the probability of short-range ordering. 

Figure 2b shows the temperature dependences of the electrical 
resistivity of polycrystalline Ago.97Cuo.03 alloy, single-crystal Ag, 
and the Ago.97Cuo.03 mixed crystal (See Figure S3 for the other sam- 
ples). The lowered resistivity at the low temperature of 10 K reflects 
the fact that the doped Cu impurities did not contribute to electron- 
impurity scattering. 

To understand the change in resistivity after the addition of Cu in 
the Ag single crystal, we first examined Matthiessen's experimental 
findings 24,25 . According to Matthiessen's rule, the electrical resistivity 
of an impure metal can be separated into two contributions, given as: 



p(c,T) = p 0 (c) + p i (T) 



(1) 



where p 0 is the temperature-independent residual resistivity, caused 
mainly by electron-impurity (el-imp) and electron-defect (el-def) 
scattering, and is dependent on the impurity concentration c, and 
pi is the temperature-dependent intrinsic resistivity, determined 
mainly by electron- phonon (el-ph) scattering 26 . For a p 0 of the binary 
solid-solution system, and considering Nordheim's rule, given as p 0 
oc c (1 - c), in which c is the mole fraction of Cu, the increased p for 
>5% Cu-doped mixed crystals can be attributed to the enlarged el- 
imp scattering. However, the change in p, must be considered in 
understanding the drop in p in the <5% Cu-doped mixed crystals. 
The intrinsic resistivity Pi can be approximated by the Bloch- 
Griineisen equation: 



T 



p(T) = a. d - ph ( — ) 



(e*-l)(l- 



x ) 



-dx 



(2) 



where 0 D is the Debye temperature, x is a variable ranging from 0 to 
0o/T, and the constant c/. d _ ph is proportional to l^Wrj/w. 2 , where X tr 
is the electron-phonon coupling constant, and co D and co p are the 
Debye and plasma frequencies, respectively 18,27 . This equation is 
favored for metal elements, such as a monovalent metals with a 
spherical Fermi surface 28 . The equation is composed of three parts: 
an electronic term (w p ), a phonon-related term(6>rj, o>d)> and an 
electron-phonon coupling term (X tt ). 

First, we consider the electronic term. The value of m p is known to 
be determined solely by the carrier concentration n and the effective 
mass of electrons m* (co p oc n/m*). Thus, the electronic contri- 
bution toward the electrical resistivity can be understood based on 
the electronic structure. The electrical resistivity (p) can be described 
by the Boltzmann transport equation, given as: 



Table 1 Optimized lattice parameters (ao), bulk modulus 

[Xtr] for Ag, Cu, and the mixed crystal 


(6), estimated Debye temperature (&d), and electron-phonon coupling constant 




ao [A] 


B [GPa] 


0d [K] 


Mr 


Ag 


4.165 


88.6 


201.8 


0.121 










0.12 31 


Cu 


3.637 


136.2 


327.1 


0.135 
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0.13 31 










0.14 32 


Cu-cu dimer. 


4.150 


90.6 


140.0 


0.109 


Cu-Cu separated 


4.150 


90.7 


203.0 


0.108 
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Figure 2 | Electrical resistivity of crystals, (a): Electrical resistivities of pure single- crystal Ag, Ag-Cu mixed crystals, and polycrystalline Ag 0 97 Cu 0 03 alloy, 
measured at 300 K, and polycrystalline Cu, single-crystal Cu, and polycrystalline Ag, measured at 293 K 19,23 . (b): Temperature dependences of electrical 
resistivity of the polycrystalline Ag 0 97 Cu 0 03 alloy, single-crystal Ag, and Ag 0 97 Cu 0 03 mixed crystal. 



p(«) = 



dfo 
de 
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where e is the charge of the carrier, N is the number of k- points used 
in the calculations, Q is the volume of a unit cell,/ 0 is the equilibrium 
Fermi-Dirac distribution function, T n /t is the relaxation time, v B> i 
denotes the group velocity, £„^ stands for the band energy, and 8 
is the delta function 29 . The subscripts n and k represent the band 
index and crystal momentum in the case of a crystal solid, respect- 
ively. v„ji can be derived from the band structure using the following 
relationship: 

1 ds„, k 
dk 



Vn,k=r- 

n 



(4) 



Thus, the value of p can be obtained from the electronic structure 
and the exact value of x. In principle, the value of x is dependent on 
the band energy e, n k and temperature; however, we assumed the 
value was constant to estimate how changing the electronic structure 
of the mixed crystal affected the electrical resistivity. Moreover, we 
assumed a rigid band approximation. The Fermi level was obtained 
from the total number of electrons in the unit cells, and was used to 
calculate p. 

The electrical resistivities were estimated from the electronic 
structures of the Ag, Cu, and mixed crystal, in which electronic 
structures of the supercell model were used for all systems. The 
model with a composition of 3.125 mol% Cu was used to estimate 
the electrical resistivity of 3 mol% Cu mixed crystals. The configura- 
tion of the Cu atoms was considered to investigate the effects of Cu- 
Cu interaction on the electronic structure and the electrical resistiv- 
ity. We simulated the two extreme systems: i) the dimerized state of 
two Cu atoms (distance between Cu-Cu (d Cu -cu) = 2.824 A), and ii) 
the furthest Cu-Cu separated system within a supercell (d Cu . Cu = 
8.300 A). The dimer state was more stable, by 12 me V, than the latter 
case (Fig. S4 in SI). The density of states of the two systems is com- 
pared in Figure S5 (in SI) and shows a clear difference between the 
systems, indicating that there is an effect of ordering on electronic 
structure. The estimated electrical resistivity from the electronic 
structure of the dimer system at room temperature can be lower, 
by —1%, than that of Ag (Fig. 3). The same value of t was used for 
all systems. These results indicated that clustering of Cu can lead to a 
reduction in the electronic resistivity, compared with pure Ag. The 



amount of reduction of p from the electronic structure was smaller 
than that of the experimental value in Figure 2. 

The lowered electrical resistivity in the mixed crystal cannot be 
fully understood by the change in the electronic structure alone. 
Thus, in addition to the electronic term, the phonon-related term, 
&d and co D , was considered to be a cause of the drop in p. From Eq. 
(2) and the equation w D = kg&o/h, where k B is Boltzmann's con- 
stant and h is Planck's constant, we can see that the Debye temper- 
ature 0 D has a significant influence on p. 0 D was calculated from the 
phonon density of states (Table 1). The 0 D value for the mixed 
crystal of the Cu-separated system was slightly greater than that of 
Ag, while that of the Cu-dimer system was much smaller. These 
results indicate that the phonon properties are also greatly depend- 
ent on Cu-ordering. 
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Figure 3 | Calculated electrical resistivity of Ag, Cu, Cu-Cu dimerized, 
and Cu-Cu separated systems. Electronic structures of the supercell were 
used to evaluate the properties. All values were obtained with the 
Boltzmann transport equation with a constant value of relaxation time of 



3.4 X 10" 



, chosen to fit the experimental resistivity of Ag. 
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It is known that ). tr is independent of the temperature and can 
be obtained from the relationship between & D , and % when the 
temperature is greater than 0.7 <9 D 33,34 , given in Eq. (5): 

h _ 2n). tr k B T 
t F th 

0.0380^ 
whereFtf,^(l ^) • 

The value of t must be known to evaluate ). tr . In this study, the 
value of t was estimated by comparing the experimental p with that 
obtained using Eq.(3). The evaluated values of A tr with t and 0 D are 
shown in Table 1 . The estimated values of Ag and Cu were consistent 
with values in the literature. The value of ). tr for the mixed crystal was 
lower than that of Ag. The reduced values were mainly due to the 
increased relaxation time, and the effect of the Debye temperature 
was limited. The relaxation time for coupling (or scattering) of the 
electrons to an acoustic phonon is dependent on the deformation 
potential, the density of materials, and the mean longitudinal sound 
velocity 35,36 . We determined that the estimated sound velocity from 
phonon dispersion and the measured density did not lead to an 
increase in the relaxation time in the mixed crystal. Thus, the reduced 
scattering (i.e., increased relaxation time) was attributed mainly to 
the reduction in the deformation potential with Cu addition. On the 
other hand, the addition of Cu decreased the density, resulting in a 
reduction in the relaxation time. Thus, as the Cu amount increased, 
the change in the deformation potential had less influence. Due to 
tradeoff between these two effects, a minimum was observed in the 
resistivity of the mixed crystal at ~3 mol% Cu. 

In addition to these intrinsic mechanisms, vacancy formation, an 
extrinsic effect, may also contribute to the low resistivity value. Our 
calculations indicated that the formation of a vacancy can be signifi- 
cantly suppressed in the presence of Cu within Ag. We estimated the 
formation energy of a vacancy though the total energy calculations of 
a cubic supercell having 32 atoms. The formation energy of a vacancy 
in Ag was estimated to be 0.74 eV, and it indicates that the concen- 
tration of vacancies can be quite high, —6.56 X 10 19 cm 3 at growth 
temperature of 900°C. An interesting finding is that when one Ag is 
replaced by Cu, the formation energy of a vacancy increases greatly, 
and is estimated to be 1.19 eV for the vacancy located at the nearest 
site of Cu and 1.13 eV for the vacancy separated from Cu. This 
indicates that a small addition of Cu into Ag can reduce the concen- 
tration of vacancies significantly. Because the vacancy should be a 
strong scattering center, this suppression may explain the reduction 
in the resistivity by Cu addition. The formation of other Cu-related 
defective states due to the low solubility limit of Cu may lead to an 
increase in resistivity, as the concentration of Cu exceeds 3 mol%. 

In summary, we observed an uncharacteristic drop in the resistiv- 
ity of single-crystal Ag grown using the Czochralski method as a 
result of adding a small amount of Cu impurity (up to 3 mol%). 
We found that the change in the electronic structure, through Cu- 
Cu dimerization, benefits the reduction in the electrical resistivity, 
resulting in the resistivity of Cu-added Ag being lower than that of 
pure Ag. However, this electronic contribution alone is not sufficient 
to explain the significant drop in resistivity. Thus, we suggest that the 
reduced scattering may be due to the reduction in the deformation 
potential with Cu addition. From an extrinsic viewpoint, the sup- 
pression of vacancy formation in the presence of Cu within Ag may 
result in a reduction in resistivity due to the reduction in the con- 
centration of vacancies, which behave as strong scattering centers. 

Methods 

Single-crystal Cu and Ag, and the Ag!_ x Cu x mixed crystals (x — 1, 2, 3, 5, and 

10 mol%) were grown using the Czochralski method. To prevent oxidation during 

crystal growth, the growth process was conducted in an Ar gas (99.999%) atmosphere. 

All samples were prepared carefully using the wire-EDM process from an ingot of 
mixed crystals 19 - 20,37 . The crystal quality of the Ag x _ x Cu x mixed crystals was investi- 
gated by pole figure measurement using high-resolution XRD (PANalytical X'pert 



PRO MRD). Radio frequency GDS (JY 10000 RF, KBSI-PA314, Korea Basic Science 
Institute, Busan Center) was used for compositional analyses of Ag x _ x Cu x mixed 
crystals and to identify any unintentional impurities. 

Electrical resistivity and temperature were measured using the four-point probe 
method and a physical property measurement system (PPMS, Quantum Design, 
USA). We kept the contact area constant using gold-coated pogo pins to minimize the 
error caused by contact resistance. We used the reversal method, in which the voltages 
for current flow in one direction and the same current flow in the reverse direction are 
measured to reduce additional voltage due to thermoelectric effects. The reverse 
method gave us reliable data by eliminating the temperature difference between the 
two readings. We used a source meter (Keithley 2425-C) as the current source and a 
nanovoltmeter (Keithley 2182A) to measure drops in voltage. 

The conventional unit cell of Ag and Cu can be described as face-centered cubic 
(FCC). For the electronic structure calculation, a primitive cell containing one atom 
was used for Ag and Cu. To describe the 3 mol% Cu-doped crystal, we constructed a 4 
X 4 X 4 supercell, and substituted two Cu atoms for Ag, representing 3. 125 mol%Cu. 
All possible configurations were considered and the difference in the total energy is 
shown in Figure S4 (in SI). The electronic structure and structural optimization were 
carried out using the Vienna ab-initio simulation package (VASP) 38 ' 39 . Projector 
augmented-wave pseudopotentials were used 40 . The exchange and correlation were 
treated within the generalized gradient approximation (PBE-GGA) 41 . The Brillouin 
zones (BZs) of the unit cell were represented by the gamma-centered Monkhorst- 
Pack scheme 42 . The BZs were sampled by at least 23 X 23 X 23 and 5X5X5 k-point 
meshes for the primitive unit cell and supercell, using the Methfessel-Paxton scheme 
with a smearing width of 0.1 eV 43 . The plane wave-cutoff energy was greater than 
400 eV for all calculations. The structures were fully relaxed using a quasi-Newton 
algorithm, until the force on each atom was within 10~ 5 eV A -1 . Electrical resistivity 
was estimated using the Boltzmann transport equation, implemented in the 
BoltzTraP code 44 , assuming a rigid band approximation and constant relaxation time. 
The values of the k-point mesh were varied to investigate the convergence of the 
electrical resistivity and an intercept value at l/(k-point mesh) — 0 from an extra- 
polation was chosen. The largest value used in the extrapolation for the supercell was 
35 X 35 X 35. 

Phonon calculations were performed using the supercell approach. Force constants 
were obtained using the framework of the density-functional perturbation theory 
(DFPT), implemented in VASP code, and phonon frequencies were calculated from 
the force constants using PHONOPY code 45 . For comparison, the phonon frequen- 
cies were also calculated using the direct Parlinski-Li-Kawazoe method, with a finite 
displacement (FD) of 0.01 A 46 . The difference between the DFPT and FD methods 
was negligible. Supercells containing 4X4X4 primitive cells (64 atoms) were used to 
calculate the phonon properties. We carried out the DFPT calculations using at least 5 
X 5 X 5 k-point grid meshes for BZ integration. 
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